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Abstract Neutrino-driven winds are thought to accompany the Kelvin-Helmholtz 
cooling phase of nascent protoneutron stars in the first seconds after 
a core-collapse supernova. These outflows are a likely candidate as 
the astrophysical site for rapid neutron-capture nucleosynthesis (the r- 
process). In this chapter we review the physics of protoneutron star 
winds and assess their potential as a site for the production of the heavy 
r-process nuclides. 

We show that spherical transonic protoneutron star winds do not 
produce robust r-process nucleosynthesis for 'canonical' neutron stars 
with gravitational masses of 1.4 Mq and coordinate radii of 10 km. We 
further speculate on and review some aspects of neutrino-driven winds 
from protoneutron stars with strong magnetic fields. 
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1. Introduction 

1.1 r-Process Nucleosynthesis 

A complete and self-consistent theory of the origin of the elements has 
been the grand program of nuclear astrophysics since the field's incep- 
tion. Of the several distinct nuclear processes which combine to produce 
the myriad of stable nuclei and isotopes we observe, none has generated 
more speculation than r-process nucleosynthesis (Wallerstcin 1997). The 
r-process, or rapid neutron-capture process, originally identified in Bur- 
bidge et al. (1957) and Cameron (1957), is a mechanism for nucleosyn- 
thesis by which seed nuclei capture neutrons on timescales much shorter 
than those for /3~ decay. The rapid interaction of neutrons with heavy, 
neutron-rich, seed nuclei allows a neutron capture-disintegration equilib- 
rium to establish itself among the isotopes of each element. The nuclear 
flow proceeds well to the neutron-rich side of the valley of /^-stability 
and for sufficient neutron-to-seed ratio (>100) the r-process generates 
the heaviest nuclei (e.g., Eu, Dy, Th, and U), forming characteristic 
abundance peaks at ^ ~ 80, 130, and 195 (e.g. Burbidge et al. 1957; 
Meyer & Brown 1997; Wallerstein et al. 1997). 

The neutron-to-seed ratio is the critical parameter in determining if 
r-process nucleosynthesis succeeds in producing nuclei up to and beyond 
the third abundance peak. If the neutron-to-seed ratio is too small, the 
r-process may only nuclei up to the first or second abundance peak. In a 
given hydrodynamical flow, the neutron-to-seed ratio is set primarily by 
the entropy (sa; 'a' here stands for 'asymptotic') of the flow, the neutron 
richness of the matter, and the dynamical timescale for expansion. The 
neutron richness is generally quantified by the the electron fraction (Y^, 
the number density of electrons per baryon). The dynamical timescale 
(xdyn) is a characteristic time for expansion, i.e. the e-folding time for 
temperature or density at a given temperature. Thus, ,Sa, Y^, and Tdyn 
set the neutron-to-sccd ratio. The higher the entropy, the lower the 
electron fraction, and the shorter the dynamical timescale, the larger the 
neutron-to-seed ratio and the higher in A the nucleosynthetic flow will 
proceed (e.g. Hoffman, Woosley, &; Qian 1997; Meyer & Brown 1997). 

1.2 Observational Motivation: A Remarkable 
Concordance 

Recent observations of neutron-capture elements in ultra-metal-poor 
([Fe/H]^ -2.5] halo stars (Sneden et al. 1996; Burris et al. 2000; McWilUam 
et al. 1995a,b; Cowan et al. 1996; Westin et al. 2000; Hill et al. 2001) 
show remarkable agreement with the scaled solar r-process abundance 
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Figure 1. The scaled solar r-process abundances (solid line) and the abundances 
from CS 22892-052 (dots with error bars) as a function of atomic number (Chris 
Sneden, private communication). 



pattern for A > 135. Particularly for atomic numbers between = 55 
and N = 75, the distribution of abundances in these halo stars is identi- 
cal with solar. Prototypical of this class are the stars CS 22892-052, BD 
-hl7°3248 and HD 115444 (Cowan & Sneden 2002). Figure 1 shows the 
scaled solar r-process abundances (solid line) together with those from 
the ultra-metal-poor halo star CS 22892-052 (points with error bars) as a 
function of atomic number. Note both the tight correspondence between 
the sun and CS 22892-052 above the second r-process peak and the large 
discrepancies at lower atomic number, in particular the elements Sn, Ag, 
and Y. The fact that a class of very old halo stars share the same rela- 
tive abundance of r-process elements above the second peak suggests a 
universal mechanism for producing these nuclei, which must act early in 
the chemical enrichment history of the galaxy. 

The fact that observations of ultra-metal-poor halo stars show sig- 
nificant and scattered deviations from the scaled solar r-process abun- 
dance pattern below A ~ 130 has implied that there are perhaps two or 
more r-process sites (Qian, Vogel, &: Wasserburg 1998; Wasserburg & 
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Qian 2000; Qian & Wasscrburg 2000). In these models of r-process en- 
richment, distinct astrophysical sites account for the observed r-process 
abundances below and above A ~ 130. Important in these consider- 
ations are the constraints set on the site or sites of the r-process by 
the integrated galactic r-process budget. For example, if all supernovae 
produce neutron stars that are accompanied by protoneutron star winds 
that produce a robust r-process, then, given a supernova rate of one ev- 
ery 50-100 years, 10~^ — 10~^ of r-process material must be injected 
into the interstellar medium per event (Qian 2000). These numbers are 
to be contrasted with those from any other potential r-process site. For 
example, if the material possibly ejected to infinity in the formation of a 
black hole from neutron-star /neutron-star mergers accounts for the total 
galactic r-process budget (Freiburghaus, Rosswog, & Thielemann 1999, 
and references therein; Rosswog et al. 1999), then lO"-*^— 10~^ of r- 
process material must be ejected per event, given an event rate of one 
every 10^ years (Kalogera et al. 2001; Qian 2000 and references therein). 

Because the total r-process budget of the galaxy is fairly well known, 
any potential site must do more than just match the entropy, electron 
fraction, and dynamical timescales required for a high neutron-to-seed 
ratio. Simply attaining the necessary physical conditions for a neutron- 
to-seed ratio above ~100 is not sufficient. The astrophysical site, given 
an event rate, must be able to consistently produce the robust r-process 
abundance pattern in ultra-metal-poor halo stars in accordance with the 
total galactic r-process budget. Thus, if M^j is the total mass of r-process 
material produced per event, the process of assessing an astrophysical 
site's potential for the r-process is merely a matter of mapping the region 
in Sa — ^e*^ — Tdyn ~ -^cj spacc physically accessible to the potential site. In 
what follows here we review the results for just-post-supernova, neutrino- 
driven protoneutron star winds. 

1.3 Protoneutron Star Winds 

The successful two-dimensional supernova explosion obtained by Bur- 
rows, Hayes, & Fryxell (1995) shows clearly a post-explosion neutrino- 
driven wind emerging into the evacuated region above the newly formed 
protoneutron star and behind the rapidly expanding supernova shock. 
Although no study was made of this outflow as a function of progeni- 
tor, such a wind phase might naturally accompany the post-explosion 
cooling epoch in many core-collapse supernovae. 

A multiple of 10^^ erg of binding energy will be lost as the nascent 
neutron star cools and contracts over its Kelvin-Helmholtz cooling time 
(~ 10 — 20 seconds) . This energy will be carried away predominantly by 
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all species of neutrino. A small fraction of that energy will be deposited 
in the surface layers of the protoneutron star, ablating material from its 
upper atmosphere and driving a hydrodynamical wind. Whether this 
wind succeeds in escaping to infinity, or is prevented by fallback and 
reverse shocks as the supernova shock encounters the density stratifica- 
tions of the overlying stellar mantle is an important future area of study. 
In addition, the actual emergence of the wind - how it overcomes the 
non-zero pressure of the matter exterior to the protoneutron star - has 
not been fully explored. In this work, however, we assume the existence 
of such outflows as evidenced by the calculation of Burrows, Hayes, &; 
Pryxell (1995) as well as those of Janka &: Miiller (1995), and Taka- 
hashi, Witti, & Janka (1994). Our goal is to explore the basic physics of 
neutrino-driven outflows and assess this site as a candidate for r-process 
nucleosynthesis. 

The basic scenario is as follows: as the supernova shock propagates 
outward in a successful explosion, the pressure in the region between 
the protoneutron star and the shock decreases and the wind, powered 
by neutrino heating, emerges into the post-shock ejecta. The surface of 
the protoneutron star is hot (temperatures of ~ 5 MeV) and has a low 
electron fraction (typically, Ye ^ 0.1). The matter there is composed 
of relativistic charged leptons, free nucleons, and trapped photons. As 
the wind is driven outward, the matter descends a gradient in density 
and temperature. The wind material is heated only within the first 
~50km, primarily via the charged-current absorption processes on free 
nucleons and i^eP ne^ , and its entropy increases con- 

comitantly. As the temperature of the matter drops below ~ 1 MeV, 
nucleons combine into alpha particles, neutrino heating ceases and the 
material expands adiabatically with entropy Sa- These charged-current 
processes also set the asymptotic electron fraction. As the matter ex- 
pands away from the protoneutron star surface and the chemical equilib- 
rium between and i>e neutrinos obtained near the neutrino decoupling 
radius (the neutrinosphere) is broken, the luminosity and energy den- 
sity of the electron-type neutrino species determine the electron fraction. 
Typically, within 10 km of the protoneutron star surface 1^ asymptotes 
to Y^. At a temperature of ~0.5 MeV (a radius in most models of 
~ 100 km) alpha particles combine with the remaining free neutrons in 
an ct-process to form seed nuclei (Woosley & Hoffman 1992). At this 
location in the temperature profile, it is the steepness of the density 
or temperature gradient that sets the dynamical timescale of the wind. 
Simply stated, short dynamical timescalcs arc favored for the r-proccss 
because for faster expansions there is less time to build seed nuclei, and, 
hence, the neutron-to-seed ratio is preferentially larger, all else being 
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equal. When the wind material finally reaches a temperature of ~0.1 
MeV (at a radius of several hundred kilometers) the r-process may be- 
gin if the neutron-to-seed ratio, as set by Y'^^ Sa, and Tdyn, is sufficiently 
high. 

Because the success of the r-process is so dependent on the neutron- 
to-seed ratio and, hence, the electron fraction, entropy, and dynamical 
timescale of the nucleosynthetic environment, we focus on these quan- 
tities and their sensitivity to various parameters of the protoneutron 
star. In addition to establishing correlations between these quantities in 
neutrino-driven wind environments, we must also fold in the constraint 
on the total mass ejected. We will show in §1.4 that this constraint on 
Mgj significantly limits the — — Tdyn space of relevance for winds. 

1.4 Previous Work 

Duncan, Shapiro, & Wasserman (1986) were the first to explore the 
physics of neutrino-driven winds. They identified a number of important 
scaling relations and the basic systematics and dependencies of the prob- 
lem. In addition, they explored the relative importance of the neutrino 
and photon luminosities in determining the resulting hydrodynamics. 
Recent investigations have focused on the potential of these outflows for 
r-process nucleosynthesis as suggested in Woosley & Hoffman (1992). 
Witti, Janka, & Takahashi (1994) and Takahashi, Witti, & Janka (1994) 
showed that although interesting a-process nuclei were created in their 
models of protoneutron star winds, conditions for a successful r-process 
fell short in entropy by a factor of ~ 5. Qian & Woosley (1996) made 
analytical estimates of the fundamental wind properties and systematics 
and compared their scalings to numbers from hydrodynamical simula- 
tions. They further explored interesting variations to their models such 
as inserting artificial heating sources and applying an external bound- 
ary pressure. These results were followed by nucleosynthetic calcula- 
tions, taking wind trajectories as the time-history of Lagrangean mass 
elements in Hoffman, Woosley, &; Qian (1997). These efforts also showed 
that with reasonable variations in the input physics protoneutron star 
winds do not achieve high enough entropy for the dynamical timescales 
and electron fractions derived. 

Cardall & Fuller (1997) extended the analytical work of Qian & Woosley 
(1996) to general-relativistic fiows. They found that significant enhance- 
ments in entropy might be obtained from compact neutron stars with 
large M/R. Several general-relativistic treatments of the problem fol- 
lowed, including the work of Otsuki et al. (2000), Sumiyoshi et al. (2000), 
Wanajo et al. (2001), and Thompson, Burrows, k, Meyer (2001). In the 
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last of these, purely transonic steady-state wind solutions were derived 
in general relativity, with a simultaneous solution for the evolution of 
the electron fraction in radius and a careful treatment of the boundary 
conditions. 

Note that we distinguish here between models of protoneutron star 
winds and bubbles. The former is typified by those works already dis- 
cussed, in which the wind velocity approaches the local speed of sound, 
or a large fraction of the speed of sound somewhere in the flow. Such 
winds were realized in the self-consistent supernova calculations in Bur- 
rows, Hayes, & Fryxell (1995). In contrast, in the work of Woosley et 
al. (1994) the protoneutron star outflow reached speeds of only a very 
small fraction of the sound speed in the region between the protoneutron 
star and the expanding supernova shock. Approximately 18 seconds af- 
ter collapse and explosion, Woosley et al. (1994) obtained entropies of 
~400 (throughout, we quote entropy in units of ke per baryon), long dy- 
namical timescales, and electron fraction (Ye) in the range 0.36 — 0.44. 
However, in their model the supernova shock reached only 50,000 km at 
these late times. In turn, this external boundary caused the wind ma- 
terial to move slowly. It remained in the heating region for an extended 
period, thus raising the entropy above what any simulation or analytical 
calculation has since obtained. Although the r-process proceeded to the 
third abundance peak in their calculation, nuclei in the mass range near 
j4 ~ 90 (particularly, *®Sr, *^Y, and ^'^Zr) were overproduced by more 
than a factor of 100. 

1.5 This Review 

In §1.2, we review the fundamental and general equations for time- 
independent energy-deposition-driven winds in Newtonian gravity and 
in general relativity. We also discuss the integrals of the flow and our nu- 
merical procedure for solving the relevant equations. In §1.3 we discuss 
some of the particulars of modeling protoneutron star winds, including 
the neutrino heating function, the equation of state, and the evolution 
of the electron fraction. We further critically examine several of the un- 
derlying assumptions that must accompany any STich model. §1.4 sum- 
marizes our results for spherical winds. In §1.5 we consider the possible 
effects of magnetic fields and speculate on a number of issues in need of 
more thorough investigation. In §1.6 we summarize and conclude. 
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2. Hydrodynamics 

2.1 The Newtonian Wind Equations 

Assuming time-independent wind solutions, the equation for mass 
conservation is simply 

V • (pv) = 0, (1) 

implying that the mass-outflow rate (M) of a wind is a constant in radius. 
In spherical symmetry, M = Airr^pv = constant. This expression yields 
a differential equation for the evolution of the matter velocity in radius, 

Idv Idp 2 
V dr p dr r 

The equation for momentum conservation, neglecting the mass of the 
wind itself, is simply 

dv IdP GM 

v— = — ^ + F^., 3 

dr p dr 

where M is the total mass of the protoncutron star. Although we in- 
clude it here for completeness, the radiation force due to the neutrinos 
(Fjy) can be safely neglected. This approximation is justified because 
the neutrino Eddington luminosity {L^'^'^ = AttGMc/ k^) is much larger 
than the neutrino luminosities that accompany the protoneutron star 
cooling/wind evolutionary phase, is the total neutrino opacity and 
is dominated by Ve,n — > pe~ and u^p ne"'" for the electron and anti- 
electron neutrino, respectively. The p- and r-neutrino opacity is dom- 
inated by neutral-current scattering off free nucleons, the wind heating 
region being unpopulated by nuclei. Including these processes, one finds 
that L^'^'^ > 10^^ erg s^^. We are thus safe in taking Fi^ = because we 
consider winds with only L^°^ < 5 x 10^^ erg s^"*^. 

Because neutrinos contribute heating and cooling to the flow, we must 
couple to these equations to the first law of thermodynamics. We define 
the net specific heating rate, q = Heating — Cooling, so that 

de . ds P dp 

dt='' = ^d-t^7Tt' 

where d/dt = [d/dt + v- V]. In the steady state, d/dt = and we obtain 



^ ds ^ dT vT dP 



fr 



We choose to reduce eqs. (2), (3), and (5) to a set of coupled differential 
equations for dv/dr, dp/dr, and dT/dr that make the physics of the wind 



Protoneutron Star Winds 



9 



sohition manifest and the solution to the problem more easily obtained. 
We start by eliminating the pressure. Expanding P differentially in p 
and T, we have that 



dP = 



dP 
dp 



dP 



6T. 



Defining (D) as 



we obtain 



^^TdP 
p dT 



(6) 



(7) 



4 + 



where Cy is the specific heat at constant volume, Cs{= dP/dp\g) is 
the adiabatic sound speed, and ct{= dP/dp\rp) is the isothermal sound 
speed. Taking eq. (6) and dropping Fj^, we can rewrite eq. (3) as 



dv 
dr 



' dP 


dp 


dP 


dT 


dp 


rp dr 


~ df 


p dr _ 



GM 



We can now eliminate dT/dr using eq. (5) so that 



dv 



dr 




dp dP_ 



T dP 
Cyp'^ df 



dp ^ q 
dr Cyv 



GM 



(9) 



(10) 



The terms dp/dr can be eliminated using eq. (2). Using eq. (8) and 
combining terms proportional to dv/dr, we obtain an expression for 
dv/dr in terms of thermodynamic quantities returned by the equation 
of state {D, P, Cy, Cg, etc.), the basic hydrodynamical variables {p, v, 
and T), and the neutrino energy deposition function (g): 



dv 
dr 




(11) 



where Ve = (2GM/r)^/^ is the escape velocity. Combining eq. (11) with 
eq. (2) and then with eq. (5), we obtain expressions for dp/dr and dT/dr, 
respectively; 



and 




(12) 



(13) 
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2.2 The General-Relativistic Wind Equations 

The time-independent hydrodynamical equations for flow in a Schwarzschild 

spacetimc can be written in the form (Thornc, Flammang, &; Zytkow 
1981; Flammang 1982; NobiU, Turolla, and Zampieri 1991) 



1 d{vy) ^ Idp 
vy dr p dr 



= 0, 



(14) 



and 



Idy 1 dP ^ 

— - -\ = 0, 

y dr £ + P dr ' 

de e + P dp q 

-r + p— = 0' 

dr p dr vy 



(15) 



(16) 



where u,.(= vy) is the radial component of the fluid four- velocity, v is 
the velocity of the matter measured by a stationary observer, 



y 



' 1 - 2GM/rc 
1 — v'^/c^ 



2\ V2 



(17) 



e (= pc? + pe) is the total mass-energy density, p is the rest-mass density, 
P is the pressure, e is the speciflc internal energy, and q is the energy 
deposition rate per unit mass. These equations assume that the mass 
of the wind is negligible, as in the Newtonian derivation. Although not 
readily apparent in the form above, eqs. (14)— (16) exhibit a critical point 
when v equals the local speed of sound. In order to make the solution 
to this system tractable and the critical point manifest we recast the 
equations in the same form as eqs. (11-13): 



dv 
dr 



v 
2t 



1 - cl/cf 
— ■u^ 



4c 



1 - v^c^ 
c1 — v"^ 



and 



D q (1- v^jc^ 



CyT y\ c1-v^ 



dp _ 2p I v"^ - vl/Ay"^ 
dr r 



p D 



— v^ 



{vy)CvTcl-v'^'' 

2 D {P + e) ( v"^ -vl/Ay'^\ 



dT 

dr rp Cv (? ' 

q ( {\ - D I (?)c\ - v'^ 



— v"^ 



Cvivy) 



A — v"^ 



(18) 
(19) 



(20) 



Protoneutron Star Winds 



11 



In the above expressions M is the protoneutron star gravitational mass 
and 

2 T dP 



D = & 



e + P dT 



(21) 
p 



Note that by taking the Hmits f/c <C 1 and Cg/c <^ 1, we recover the 
Newtonian wind equations in critical form. 

2.3 Conservation and Numerics 

There are two integrals of the flow which we use to gauge the accuracy 
of our solution to the Newtonian or general-rclativistic wind equations. 
The first is the mass outflow rate, obtained from direct integration of the 
continuity equation (eq. 14), which yields the eigenvalue of the steady- 
state wind problem, M = Airr^pvy (y = 1 in the Newtonian limit). The 
second is the Bernoulli integral, modified by energy deposition. In the 
Newtonian case, 

MA e + -t;2 + --^ = / <fr'pq = Qir), (22) 

where Ry is the coordinate radius of the protoneutron star surface and 
the A expresses the change in the quantity in parentheses between Ry 
and r. We retain the subscript on to emphasize that throughout 
this work the coordinate radius of the protoneutron star is assumed 
to coincide with the radius of decoupling for all neutrino species. In 
general relativity, with g = 0, 'yh^/—gQo is a constant. Here, 7 is the 
Lorenz factor and h is the specific enthalpy. With a source term, the 
differential change in neutrino luminosity is given by 

e-^^|^(L.e^'^) = -q, (23) 

where dji/dr = AirPpe^. The expression ds'^ = —e^'^dt'^ + e^^dr'^ + r'^dQ 
defines the metric. The total energy deposition rate is then, 

roo 

Q = A'K drr'^ pqe^e^^. (24) 



We solve the system of wind equations in critical form using a two- 
point relaxation algorithm between Ri, and the sonic point on an adap- 
tive radial mesh (see Thompson, Burrows, & Meyer 2001; Press et 
al. 1992; London &; Flannery 1982). We employ physical boundary con- 
ditions for the transonic wind problem, enforcing thermal (g = 0) and 
chemical equilibrium (dYg/dr = 0; see §1.3) at R^, and u = at the 
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outer edge of the computational domain. The last required boundary 
condition, which closes the system of equations, is an integral condition 
on the electron-neutrino optical depth (r): Ti,^{R^) =2/3. This condi- 
tion is combined with the other two at Ri, in a triple Newton-Raphson 
algorithm so that all boundary conditions are satisfied simultaneously. 
For modest radial zoning we typically maintain constant M and consis- 
tent Bernoulli integral to better than 1%, having imposed neither as a 
mathematical constraint on the system of wind equations. 

3. Particulars of Protoneutron Star Winds 

3.1 Electron Fraction Evolution 

The charged-current electron- type neutrino interactions on free nucle- 
ons - fgn <^ e~p and UgP e+n - affect the evolution of the electron 
fraction in radius; 

(vy) =Xn [T^^n + Te+n] - [^P,p + ^e-p] , (25) 

where X„ and Xp are the neutron and proton fraction, respectively. 
The r's are the number rates for emission and absorption, taken from 
the approximations of Qian Sz Woosley (1996). The subscripts denote 
initial-state particles. Because the resulting nucleosynthesis in a given 
wind model is so sensitive to the Y^. at the start of the r-process, we 
solve the coupled system of the wind equations with dYe/dr simultane- 
ously. Ignoring the details of transport and neutrino decoupling near 
the neutrinosphcrcs, the asymptotic electron fraction Y^ is determined 
by both the luminosity ratio Lp^/L^,^ and the energy ratio {spe) 
where (s,^) = {El)/{Ey), and En is the neutrino energy. To rough ap- 
proximation (Qian et al. 1993; Qian &; Woosley 1996), 

..a r,,n ^ / L,^ (gPe) - 2A + 1.2AV(epJ 

' r,,„ + rp,p \^ L,, (£,J + 2A + 1.2AV(£,J 

where A(= rUn — rup 2± 1.293 MeV) is the energy threshold for the Ue 
neutrino absorption process, v'eP ~* ne^. There are several important 
effects in protoneutron star winds which act to increase Y^: (1) The 
Threshold Effect: even if is constant in time, if both {spj 

and (ffi/g) decrease, Y^ must increase as a result of A and (2) The Alpha 
Effect: as the flow cools in moving away from the protoneutron star and 
a-particles are preferentially formed residual excess neutrons will capture 
Uf, neutrinos, thus increasing Y^ (Fuller Sz Meyer 1995; McLaughlin, 
Fuller, & Wilson 1996). The a-effect is particularly important for flows 
with long dynamical timescales. 
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3.2 Neutrino Energy Deposition 

The Charged Current Processes: In the protoneutron star wind 

context, the charged-current processes (z/gn ^ e~p and P^p ^ e+n) 
compete with neutrino-electron/positron scattering as the dominant en- 
ergy deposition mechanisms. Ignoring final-state blocking and assuming 
relativistic electrons and positrons, the charged-current specific cooling 
rate can be written as 



Ccc ^ 2.0 X 10^^ 



(27) 



where F„(y) = x"(e'^'"2' + 1)"^ dx, T is in MeV, and rje = He/T. The 
heating rate is 



i^ec - 9.3 X WR-i 



Xr^Llliel) + X^Lfjel) <^<^E{r), (28) 



where R,ye is the neutrinosphere radius in units of 10^ cm, Ly{= 10^^ L^^) 
and (4) are defined at R^, $ = [(1 - 2GM/R^c^)/{1 - 2GM/rc^)Yl'^ 
is the gravitational redshift, and E!(r) is the spherical dilution function. 
In the vacuum approximation and assuming a sharp neutrinosphere, 

S(r) = 1 - ^1 - (i?./r)V$2. (29) 

The redshift term, appearing in eq. (29), accounts for the amplifi- 
cation of the heating processes due to the bending of null geodesies in 
general relativity (Salmonson & Wilson 1999). Cardall & Fuller (1997) 
showed that although this amplification is important, the dominant 
general-relativistic effect on the heating rates is due to the gravitational 
redshift of the neutrino energy and luminosity (note the in cq. 28) . 

With eqs. (27) and (28), the specific energy deposition rate (ergg~^ s~^) 
due to the charged-current processes is simply = Hcc — Ccc- 

Inelastic Neutrino Scattering: At high entropies, electron-positron 
pairs are produced in abundance. The energy transfer associated with a 
single neutrino-electron or positron scattering event is well approximated 
by u>i ~ {e^^ — 4T)/2, where u> is the energy transfer and i labels the 
neutrino species (Bahcall 1964). Note that to allows for both net heating 
or net cooling, depending upon the local temperature and the neutrino 
spectral characteristics. A number of researchers have dropped the '4T' 
coofing part from uj (Qian & Woosley 1996; Otsuki et al. 2000). This 
leads to larger net energy deposition and significant modifications to the 
wind solution. We have made a comparison between models with and 
without this term and find that omitting cooling leads to a 40 — 60% 
increase in M and a 10 — 25% decrease in Tdyn- The asymptotic me- 
chanical luminosity, Mv'^/2, is increased by as much as 80% and the 
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asymptotic entropy is decreased at the 5% level. The largest deviations 
are in the models with the lowest luminosity and the highest Sa- This 
is expected because as the luminosity decreases and the entropy of the 
flow increases, neutrino-electron scattering contributes more to the total 
energy deposition profile. Because the shape of the overall heating pro- 
file is not significantly affected by omitting the '4r' cooling term from 
cj, the asymptotic entropy is only modified slightly. 

For inelastic neutrino-electron scattering, the net specific energy de- 
position rate can be approximated by g ~ cn(>ny.{ay.eOo), where rie and 
n^.. are the number density of electrons and neutrinos, respectively, and 
ay^e c± KiTe^i (Tubbs and Schramm 1975), Ki = ao^i/2ml is a neutrino 
species dependent constant, where rrif, is the mass of the electron in MeV, 
(Jo ~ 1.71 X 10~^^ cm^, and Aj is the appropriate combination of vector 
and axial-vector coupling constants for neutrino species i. Averaging 
properly, we find that 




F2M 



$2 gj.g g 1 g 1 



AT 



FaM 



(30) 



where rji, is an effective neutrino degeneracy parameter (Janka and Hille- 
brandt 1989), and is the flux factor, which is related to eq. (29) by 
(fi) = /?^/2<l>^H(r)r^. In order to obtain the contribution to the net 
heating from neutrino-positron scattering, rje —rje and one must also 
make appropriate changes to Aj. 

Electron/Positron Annihilation and its Inverse: Also at high 
entropies, cooling and heating due to e+e~ ^ fjMj must be included. As- 
suming relativistic electrons and positrons, and ignoring Pauli blocking 
in the final state, the cooling rate is 



C~ 1.4 X IO^^tVs V(%) ergsg-^s-i, 



(31) 



where f{r]e) = [F4(r/e)i^3(-%) + i^4(-??e)i^3(??e)]/2F4(0)F3(0), P8 is the 
mass density in units of 10^ g cm~^ and T is in MeV. The specific heating 
rate due to the inverse process, i^P e+e~, is simply (Qian & Woosley 
1996) 



i7 ~ 1.6 X 10^^ 
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, (32) 



where *(x) = (1 - x^ix"^ + Ax + 5), and x = (1 - (i?^/r)V$2-)i/2_ 
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3.3 Equation of State 

At the temperatures and densities encountered in protoneutron star 
winds, exterior to the radius of neutrino decouphng, to good approxi- 
mation, free neutrons, protons, and alpha particles may be treated as 
non-relativistic ideal gases. We also include photons and a fully general 
electron/positron equation of state is employed (Evonne Marietta, pri- 
vate communication). Sumiyoshi ct al. (2000) have found that using a 
general electron/positron EOS can decrease the dynamical timescale in 
the nucleosynthetic region of the wind by as much as a factor of two. 
Such a modification is important when considering the viability of the 
neutrino-driven wind as a candidate site for the r-process. For this rea- 
son, a general electron/positron EOS is essential. Also of importance 
is the inclusion of alpha particles. The formation of alpha particles ef- 
fectively terminates energy deposition via the processes VgU •s-^ e~p and 
i>eP ^ e^n. Failure to include alpha particles results in more heat- 
ing and a broader energy deposition profile. The entropy of the flow is 
thereby higher. In low luminosity, late-time, high entropy (sa ~ 200) 
the difference in Sa is ~ 20 units. 

4. Results: Spherical Models 

For a given M, and Lj,, the solution to eqs. (18)-(20) yields radial 

profiles of temperature, density, electron fraction, and velocity. From 
these quantities, one obtains Sa, Y^, and Tdyn - the critical parameters 
in determining the neutron-to-seed ratio and, hence, the resulting r- 
process. By comparing these numbers with r-process nucleosynthesis 
survey calculations like those of Hoffman, Wooslcy, & Qian (1997) or 
Meyer &: Brown (1997), one can quickly see if a given wind solution 
inhabits a point in Sa-^e^-^dyn space where a robust 3'''^-peak r-process 
is likely. 

In order to map the wind solution space, we have constructed evo- 
lutionary tracks from our steady-state models. As the supernova com- 
mences, we expect the protoneutron star to have large radius and high 
neutrino luminosity. As the cooling epoch proceeds, the protoneutron 
star will contract to its final radius (perhaps 10 km) and the luminosity 
may decrease as a power-law or quasi-exponentially in time (Burrows k. 
Lattimcr 1986; Pons et al. 1999). The actual time dependence of Ly and 
Ri, depends sensitively on the equation of state of dense nuclear matter 
and the details of transport and deleptonization of the protoneutron star 
by neutrinos. 

Figure 2 summarizes the results of Thompson, Burrows, &: Meyer 
(2001) for protoneutron stars with gravitational masses of 1.4, 1.6, 1.8, 
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and 2.0 Mq. Each thin sohd line is a sequence of steady-state gcneral- 
relativistic wind models in the plane of Tdyn versus Sa- Each track, for a 
given M, starts with R ~ 20 km and = 8.0, corresponding to a total 
neutrino luminosity of 3.7 x 10^^, where L^^ = L,^/10^^. The models 
with largest Ljy and Ry have the lowest Sa (~ 50 — 70) and moderate 
'^dyn (~ 9 vas). We take Ljy oc t~^''^ and Rv{t) such that the protoneutron 
star radius decreases linearly in time from 20 km to 10 km in one second 
(thin solid lines labeled, 'Fast Contraction'). For comparison, we also 
include a model that has Ry{t) oc t^^/^ for M = 1.4 (labeled, 'Slow 
Contraction'). In the 'Fast' cases, all models move to much higher Sa 
as Rv goes from 20 km to 10 km. The reached at each Ly is set in 
part by M/Ry{t), with the 2.0 model reaching Sa — 150 when Ry 
reaches 10 km. Once the protoneutron star reaches 10 km, R^ is fixed 
and each track makes a sharp turn toward much longer Tdyn and only 
moderately higher Sg,. Due to the relatively slow contraction, the model 
with Ru{t) oc t~^^^ never exhibits such a sharp turn in the Sa~7"dyn plane 
and eventually joins the 'Fast' 1.4 evolutionary track at Tdyn ~ 0.015 
seconds. Note that after R^ has reached 10 km, the evolutionary tracks 
evolve along characteristic curves in the Sa-Tdyn plane. We find that 
these two r-process parameters approximately follow the power law 

oc T°^„. (33) 

In their analytic and Newtonian exploration of protoneutron star winds, 

1 /6 

Qian & Woosley (1996) found that Sa T^y^^ at constant R^ and M. 
This 1 / 6*^^ scaling and that found in eq. (33) , from our general-relativistic 
wind solutions, are to be compared with the analytic work of Hoffman, 
Woosley, & Qian (1997). They find (sec their eqs. 20a & 20b) that the 
Sa required to achieve 3'*^-peak r-process nucleosynthesis as a function 
of Tdyn, at constant Y^, is given by 

Sa ~ 730 ( ^) (^^) kB baryon-i (34) 

vo.soy vo.io s; ^ ' 

for > 0.38. In Fig. 2 we show (thick solid lines) results from eq. (34) 
for Y^ = 0.48 and 0.38. Because all of the several hundred wind mod- 
els in Fig. 2 have 0.45 ^ Y,^ ^ 0.495, we conclude that models with 
M = 1.4 Mq fall short of the required entropy by at least a factor of 
~ 3. These lines are meant only to delineate the relevant range of Sa 
and Tdyn required for 3'^'^-peak nucleosynthesis. Actual nucleosynthetic 
calculations in the wind profiles themselves are preferred to the sim- 
ple comparison on this plot. Thompson, Burrows, & Meyer (2001) did 
just this, carrying out the full r-process calculation in the 1.4 M© evolu- 
tionary track. They found that nucleosynthesis did not proceed beyond 
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Figure 2. Evolutionary tracks of steady-state wind models showing the various 
correlations between Sa and Tdyn- The solid lines show models for 1.4, 1.6, 1.8, and 
2.0 Mq protoneutron stars (from Thompson, Burrows, & Meyer 2001). The thick 
solid linos show an analytic approximation from Hoffman, Woosley, & Qian (1997) 
for the Sa required for 3'^'^-pcak r-proccss nucleosynthesis for = 0.48 and 0.38. If 
Sa goes above this line then production of the S'^'^-peak nuclides is likely. The thick 
dashed line denotes the rdyn beyond which, for given M, M is too small for the wind 
to contribute significantly to the total galactic r-process budget. 



A ~ 100, in accordance with the predictions of the survey calculations 
of Meyer & Brown (1997) and Hoffman, Woosley, & Qian (1997). 

That the exponent 1/3 appears in eq. (34) and 0.2 appears in eq. (33) 
implies that once Ri, is set, for a given M, the wind cannot evolve 
into a region of the Sa-Tdyn plane where the r-process can take place. 
Put another way, the nearly horizontal lines of constant Ry, for a given 
M, in Fig. 2 cannot cross the lines of cq. (34) at arbitrarily long Tdyn- 
Thus, we conclude, along with Takahashi, Witti, & Janka (1994), Qian 
& Woosley (1996), Suniiyoshi et al. (2000), Wanajo et al. (2000), and 
Otsuki et al. (2000), that winds from protoneutron stars with 'canonical' 
parameters M = 1.4 M© and Ri, = 10 km fail to produce robust r-process 
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nucleosynthesis up to and beyond the 3'' r-proccss peak. Although the 
effects of general relativity arc important in determining the dynamical 
timescales and entropy of the wind as predicted in Cardall &: Fuller 
(1997), for reasonable Y^'s, the entropy falls short of that required for 
S'^^-peak nucleosynthesis by a large factor (~ 3, slightly better than the 
factor of ~ 5 found by Takahashi, Witti, & Janka 1994). 

4.1 Mass Loss 

Having solved for M at every point along these evolutionary tracks, 
and assuming that L^{t) oc t~^'^, we can calculate the total mass ejected 
as a function of time: 

Mej(t) = /* M{t')dt'. (35) 

Observations of ultra- metal-poor halo stars (see §1.1.2) suggest that the 
astrophysical site for production of the heavy r-process nuclides is uni- 
versal and acts early in the chemical enrichment history of the galaxy. 
The fact that if all supernovae produce this r-process signature, then 
10~^ — 10~^ M© of r-process material must be ejected per event (Qian 
2000) allows us to constrain the space of relevant wind solutions with 
eq. (35). The heavy dashed line in Fig. (2), assuming Ly{t) oc t~^'^ , 
shows the point along each evolutionary track in the Sa-Tdyn plane be- 
yond which (to longer Tdyn) M is simply too small to generate the re- 
quired total mass loss so as to contribute significantly to the total galac- 
tic r-process budget. If r-processing begins to the right of this line, 
less than 10~^ M© will be ejected. Although the position of this AMgj 
line must change for different L,j{t) and Ru{t), such a bound must ex- 
ist for any cooling model. Thus, although the wind may eventually 
evolve to arbitrarily long dynamical timescales, we conclude that the 
range of Tdyn relevant for r-process nucleosynthesis is significantly con- 
strained by consideration of M and Mgj. For example, the track for 
M = 1.4 M0 reaches Sa ~ 400 only when Tdyn is several seconds and M 
is of order 10~^^ Mq s~^. Hence, even if the r-process could commence 
in this epoch, it would need to persist for ~ 10^ seconds in order to 
contribute significantly to the total galactic r-process budget. Conser- 
vatively, then, if transonic protoneutron star winds are the primary site 
for the r-process, this constraint on the amount of mass ejected per su- 
pernova implies that the epoch of r-process nucleosynthesis must occur 
for Tdyn less than ^ 0.07 — 0.1 seconds. 
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4.2 Conclusions from Spherical Models 

Models of transonic winds from neutron stars with M ~ 1.4 Mq and 
Ru ~ 10 km fail to produce the heavy r-process nuclides. Models with 
much higher gravitational mass and even smaller coordinate radii, with 
large neutrino luminosities can achieve 3''^-peak nucleosynthesis. We 
find that with M ~ 2.0 M© and R^, ~ 9 km that some S'^'^-peak nuclides 
are produced. It is difficult to understand how such massive and com- 
pact objects might be created in standard supernova scenarios. This 
has led Thompson, Burrows, & Meyer (2001) to speculate that the near 
environments of collapsars, black holes caused by stellar collapse, sur- 
rounded by a thick accretion disk, might generate outflows much like 
the wind solutions described here, but benefiting from the general rela- 
tivistic effects as with a M ~ 2.0 M© and Ri, ~ 9 km protoneutron star. 
Barring these possibilities, however, we are left with possible modifica- 
tions to the physics described here which might lower Y^, increase Sa, 
or decrease Tdyn- 

Some possibilities for decreasing include (1) different i^e and 
spectral characteristics (see eq. 26), (2) neutrino oscillations (Qian & 
Fuller 1995a,b), and (3) neutrino transport effects. Of relevance for (1), 
changes to the high-density nuclear equation of state may interestingly 
effect the electron and anti-electron neutrino spectra. Modifications to 
the neutrino energy-deposition profile may effect both Sa and Tdyn- An 
increase in q{r) at fairly large radius (50 — 100 km) can increase Sa and 
decrease Tdyn (Qian & Woosley 1996). Such a modification to q{r) could 
be caused by non-standard neutrino physics or even by magnetic field 
reconnection (see §1.5). 

Finally, one may suspect that the assumption of sphericity is a funda- 
mental problem - that protoneutron star winds and their ejecta simply 
cannot be understood fully in one spatial dimension. In particular, one 
may add additional degrees of freedom and break spherical symmetry, 
by considering the effects of rotation and magnetic fields. As a start 
to the very complex problem of full magnetohydrodynamic and gencral- 
relativistic outfiows, in the following we attempt to quantify some of the 
basic numbers and scalings. 

5. Magnetic Protoneutron Star Winds 

The solar wind cannot be explained in detail without consideration 
of magnetic effects on the outflow. This, coupled with the fact that a 
class of neutron stars are observed to have very high surface magnetic 
fleld strengths {magnetars, Kouveliotou et al. 1999; Duncan &; Thomp- 
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son 1992), motivates an examination of MHD effects on neutrino-driven 
protoneutron star winds and their nucleosynthetic ejecta. 

To date, such effects have received little attention. 

Qian k. Woosley (1996) speculated quahtatively on the role of mag- 
netic fields in their wind solutions, noting that tangled field topologies 
might impede the flow in escaping to infinity, but the effects discussed 
were not quantified. Nagataki & Kohri (2001) considered a monopole- 
like magnetic field with rotation in one dimension by restricting their 
attention to the equatorial plane. The formulation was directly analo- 
gous to that of Weber &; Davis (1967). However, because of the complex 
critical point topology encountered in this MHD wind problem, they 
were unable to assess the importance of field strengths above ~ 10^^ G, 
although they were able to consider a variety neutron star rotation peri- 
ods. Recently, Cameron (2001) argued qualitatively that core collapse, 
rotation, and magnetic fields conspire to form jets and a post-collapse 
accretion disk that feeds these jet outflows. 

Below we briefly discuss a small subset of the possible effects expected 
from neutrino-driven MHD winds. 

5.1 Non-Spherical Expansion 

In a strong magnetic field, the character of the neutrino-driven outfiow 
may be significantly modified by the non-spherical divergence of open 
field lines, along which the wind is channeled. Kopp k. Holzer (1976) 
first considered these effects in their models of the solar wind in coronal 
hole regions. In this case, eq. (1) becomes 

±{Apv) = 0, (36) 

where ds is the differential line element along the magnetic field and A{s) 
is an arbitrary area function. The derivation of equations analogous to 
eqs. (11)-(13), starting with eq. (36), is straightforward. The solution 
to those equations, however, is complicated by the possibility that more 
than one critical point may exist in the flow if A{^s) changes rapidly 
(Kopp k, Holzer 1976; Bailyn, Rosner, & Tsinganos 1985). Indeed, 
standing shocks, connecting physical solutions, may exist in the flow 
for rapid areal divergence (Habbal k. Tsinganos 1983; Bailyn, Rosner, k 
Tsinganos 1985; Leer k Holzer 1990). For smoothly and modestly chang- 
ing A{s), however, the solution proceeds as in the spherical case - the one 
dimensional problem now along ds instead of dr. Charboneau k Hund- 
hausen (1996) have constructed quasi-two-dimensional models of flow in 
the field lines of the open region in a helmet streamer configuration (see 
also Pneuman k Kopp 1971; Low 1986). Helmet streamer /coronal hole 
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magnetic field configurations in the context of the sun have been the 
focus of considerable theoretical effort (e.g. Mestel 1968; Pncunian & 
Kopp 1970; Pneuman & Kopp 1971; Steinolfson, Suess, & Wu 1982; Us- 
manov et al. 2000; Lionello et al. 2002). In these models, pressure forces, 
inertia, gravity, and a strong ordered dipole magnetic field conspire to 
produce a region of closed magnetic field lines close to the central star, at 
latitudes near the magnetic equator. At the magnetic poles, the flow is 
radial. At intermediate latitudes, between the pole and the closed zone, 
open magnetic field lines bend toward the equator close to the star and 
then extend radially. In these models A(s) exhibits smooth variations 
and only for streamlines emerging from latitudes very near the closed 
zone do large deviations from purely radial flow exist (Charboneau &; 
Hundhausen 1996). 

Taking yl(r) = r'^f{r) and/(r) = (/max exp[(r-iii)/o-]+/i)/(exp[(r- 
Ri)/a] + 1), where /i = 1 - (/^ax - l)exp[{R„ - Ri)/a], as in Kopp 
& Holzer (1976), we have computed several Newtonian wind models for 
comparison with purely spherical expansion. This function varies most 
rapidly near with the change in /(r) over a radial distance i?i ± a. 
With /max = 4 so that A{r) is four times as large at a given radius as 
spherical expansion and with Ri = 2Ri, and a = Ri, so that the diver- 
gence is smooth, we find that Tdyn increases by ~ 60%, M decreases by 
~ 65%, Sa increases from 68 to 74kB baryon"^, and that the asymptotic 
mechanical luminosity (Pmcch = ^^^/2) drops by more than a factor of 
three. In contrast, constricting the flow with /max = 1/4 yields a much 
faster wind. In this C&iSG, ^dyn 

decreases from 3.2 milliseconds to 1.7 
milliseconds. Sa decreases from 68 to 66.5 ke baryon^^. M and Pmcch 
increase by a factor of ~ 3 and ~ 5.5, respectively. Clearly the quality 
of the areal divergence can significantly infiuence the properties of the 
flow. Although Sa was not affected by more than a few percent, the 
changes in Tdyn evidenced by this simple comparison imply that a more 
thorough investigation is warranted. We save a detailed exploration of 
these effects for a future work. 

5.2 Closed Zones &; Trapping 

The ideas of this section have recently been set down in Thompson 
(2003). Here, we follow the discussion of Thompson (2003) closely. 

Figure 3 shows profiles of temperature, entropy, energy deposition 
rate, and pressure for a high neutrino luminosity protoneutron star wind 
model with M = 1.4 Mq and = 10 km. Also shown (thick dot- 
ted line) is the magnetic energy density i?^/87r. For simplicity we take 
B = Bq^Rb/t)^-, where Rb is a reference radius for the magnetic field 
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footpoints. Because of the exponential near-hydrostatie atmosphere in 
these wind models (note sharp drop in P in Fig. 3), we take Rb = 11 km. 
The surface magnetic field strength Bq is here set to 1.5 x 10^^ G. We 
define the quantity (3 = P/{B^/8-7t) and Rp as the radius where /3 = 1. 
Here, R/^ ~ 46.5 km. From this figure it is clear that a ~ 10^^ G field 
can dominate the matter pressure during the wind epoch. Note that this 
figure is only for a single wind model, with a single neutrino luminosity. 
As Li, drops, the pressure profile drops everywhere so that R^ moves out 
in radius for constant Bq. Because P{r) drops everywhere as L^, decays, 
at any instant in time, a lower Bq is required such that /3 ~ 1 at some ra- 
dius. This implies that even though a 10^^ G field may not dominate the 
wind dynamics at early times, if the protoneutron star somehow main- 
tains this field strength as drops, the field will eventually dominate 
as P{r) decreases (Thompson 2003). 

Prom this admittedly limited comparison, we conclude that neutrino- 
driven winds from protoneutron stars with magnetar-like surface field 
strengths may be significantly affected by the presence of such a field. 

Prom MHD models and observations of the solar wind, we expect that 
a strong magnetic field, that dominates the wind pressure inside Rp, may 
form a closed zone at these radii, at latitudes near the magnetic equator. 
The configuration of the flow would then be analogous to the helmet 
streamer described in §1.5.1 (Steinolfson, Suess, & Wu 1982; Usmanov 
et al. 2000; Lionello et al. 2002). If heating and cooling balance so 
that q = throughout the closed zone, this structure may be stable 
and the matter in this regime permanently trapped. However, if heating 
dominates then the pressure of the trapped matter will increase. Thus, if 
the material is trapped inside Rfs with P <^ i^^/Bvr, net neutrino energy 
deposition in the closed zone (see Fig. 3) must increase P to ~ B^/8tt. 
If this happens, we expect the matter to escape dynamically. In this 
way, the closed magnetic field structures that form where f3 < 1 are 
unstable. Importantly for the r-process, because the pressure increase 
of the matter is caused by neutrino heating, it is necessarily accompanied 
by an increase in the matter entropy (Thompson 2003). 

Very roughly, the matter will be trapped for a time set by P, B'^/Stt, 
and q; 



Assuming that T, p, and q do not change significantly in Ttrap, the en- 
tropy amplification associated with such an increase in pressure is then 
(Thompson 2003) 



Ttrap - [5V87r-P]/[gp]. 



(37) 



(38) 
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Figure 3. Temperature T (McV, dotted line), entropy s (10 ke baryon ^, thin solid 
line), energy deposition rate q (lO^'^ erggs"^, long dashed line), log of the pressure 
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^, thick solid line), and log of the magnetic energy density S^/Stt 

short dashed line) for a protoneutron star wind model with M = 
1.4 M© and Ri, = 10 km. Here, we take B = Bo{RB/r)^ , where Rb = 11 km and 
Bo = 1.5 X IQi^ G. 



When eqs. (37) and (38) are evaluated at a characteristic radius for en- 
ergy deposition (say the half-asymptotic-entropy radius), they yield an 
order of magnitude estimate for As. For very high Bq and low P{r) (slow 
winds with low Li,), Thompson (2003) found that there is a radius (Rq) 
inside of which cooling balances heating before P approaches B'^/Stt. 
Therefore, in this simple picture, the matter interior to Rq is perma- 
nently trapped (barring MHD instabilities that might very well arise). 
Importantly, for any Bq, Rq is always less than Rp so that the trapped 
matter between these two radii can escape with As set by eq. (38). 

Figure 4 shows Sa versus Tdyn (analogous to Fig. 2) for a large set of 
protoneutron star wind models (Thompson 2003). The thick solid line 
shows spherical, steady-state models as described in §1.4 for constant 
= 10 km and M = 1.4 M©. The dashed line shows again the analyt- 
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Figure. 4- Non-magnetic, spherically symmetric wind models (thick solid line) in 
the plane of Sa versus Tdyn for protoneutron star wind models with M = 1.4 and 
Ri, = 10 km, for a large range in neutrino luminosities. Thin solid lines show Sa as 
a function of Tdyn, employing the entropy enhancement as in eq. (38) for 2 x 10^'' < 



Bo < 2 X lO^'^ G. The dashed hne shows the results of eq. (34) for 
Hoffman, Woosley, & Qian 1997). Figure from Thompson (2003). 



0.48 (from 



ical results of Hoffman, Woosley, 8z Qian (1997) for the required, at 

a given T^yn, for 3'''^-pcak r-proccss nucleosynthesis (see eq. 34; compare 
with Fig. 2). Above this dashed line, for = 0.48, the neutron-to- 
seed ratio is high enough for a robust r-process. The thin solid lines 
are constructed from the non-magnetic wind models (thick solid line), 
by applying eqs. (37) and (38) for Bq = 2 x 10l^ 2 x 10l^ 4 x 10^^, 
6 X 10^"^, 8 X 10^^, 1 X 10^^ and 2 x 10^^ G. From Fig. 4 it is clear both 
that the spherical non-magnetic wind models fall short of the entropy 
required for 3'^'^-peak r-process and that for Bq ~ 10^'^ G, the entropy 
enhancements caused by trapping may be sufficient to account for this 
deficit. 

There arc a number of effects that might decrease the entropy en- 
hancements discussed here as a result of trapping in closed magnetic 
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field structures. Any physical effect that globally disrupts the closed 
zone on a timescale much less than rtrap would significantly undermine 
the entropy enhancements estimated in this scenario. Such effects might 
include MHD instabilities, differential rotation, and rapid motion of the 
magnetic field footpoints due to convection (Thompson 2003; see also 
Duncan & Thompson 1992; Thompson &; Duncan 1993; Thompson &; 
Murray 2001). It is worth noting that the very early configuration of 
the protoneutron star magnetic field is highly uncertain and may be a 
complex of high-order multipoles. The large-scale closed zone described 
in Thompson (2003) might not then obtain. However, in this case many 
closed zones may exist and eq. (38) may be used to estimate the entropy 
enhancement locally at many sites on the surface of the protoneutron 
star as closed regions with a variety of /3 are generated and then opened 
by neutrino heating. If the field is very complex, twisted, or sheared 
reconnection may deposit energy in the flow, as first suggested by Qian 
& Woosley (1996) in this context. Studies of extra energy deposition 
show that this may increase or decrease the entropy of the flow, depend- 
ing crucially on where the energy is deposited (Qian & Woosley 1996; 
Thompson, Burrows, Sz Meyer 2001). 

6. Summary, Conclusions, Sc Implications 

The subject of protoneutron star winds is relatively young, born only 
in the early-1990s. Much of the physics attending the emergence and 
evolution of these outflows in the just-post-supernova environment is un- 
certain and intimately tied with other outstanding issues in neutron star 
birth: rotation and magnetic fields. Many of the phenomena intensively 
investigated in the context of the sun reconnection, flares, coronal 
mass ejections, closed loops, prominences, flux emergence, coronal holes, 
streamers - may play important roles in determining the nucleosynthetic 
consequences of the wind/cooling epoch. 

The results presented in §1.4 from Thompson, Burrows, & Meyer 
(2001) as well as the results of Takahashi, Witti, & Janka (1994), Qian 
& Woosley (1996), Sumiyoshi et al. (2000), Otsuki et al. (2000), and 
Wanajo et al. (2000) indicate that spherical steady-state winds from 
canonical neutron stars cannot attain the requisite entropy for robust 
r-process nucleosynthesis. Figure 2 shows, however, that very com- 
pact, massive, and luminous neutron stars may realize a short dynamical 
timescale, modest entropy r-process. In addition, it may also be that the 
actual electron fraction of protoneutron star winds is much lower than 
that derived in Thompson, Burrows, Sz Meyer (2001) and implied by 
the neutron star cooling calculations of Pons et al. (1999). If could 



26 



be made to be ~ 0.3 in the models of Fig. 2, the 1.4 Mq evolutionary 
track might naturally generate 3'^'^-peak nuclides. Then again, perhaps 
1.4 Mq neutron stars are responsible for production of only the r-process 
elements below A ~ 130 as in Fig. (1) (see §1.1.2). 

Section 1.5 shows that magnetic effects can change the entropy and 
dynamical timescale of a given flow solution considerably for surface field 
strengths of order 10^^ G. It may be that only neutron stars born with 
magnetar-like field strengths produce robust r-proccss signatures. In any 
case, the need for multi-dimensional magnctohydro dynamic simulations 
of wind emergence and evolution are required to address many of these 
still open questions more fully. 
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